arXiv:astro-ph/9706159v2 16 Dec 1997 


A New Approach to Probing Large Scale Power with Peculiar 

Velocities 

Hume A. Feldmanj^ 

Department of Physics & Astronomy, University of Kansas 
Lawrence, KS 66045 

Richard Watkins^ 

Department of Physics & Astronomy, Dartmouth College 
Hanover, NH 03755 

Received_; accepted_ 


^feldman@ukans.edu 
^watkins@mojo.dartmouth.edu 





2 


ABSTRACT 


We propose a new strategy to probe the power spectrum on large scales 
using galaxy peculiar velocities. We explore the properties of surveys that cover 
only two small helds in opposing directions on the sky. Surveys of this type 
have several advantages over those that attempt to cover the entire sky; in 
particular, by concentrating galaxies in narrow cones these surveys are able to 
achieve the density needed to measure several moments of the velocity held 
with only a modest number of objects, even for surveys designed to probe scales 
> 100/i“^Mpc. We construct mock surveys with this geometry and analyze 
them in terms of the three moments to which they are most sensitive. We 
calculate window functions for these moments and construct a statistic which 
can be used to put constraints on the power spectrum. In order to explore the 
sensitivity of these surveys, we calculate the expectation values of the moments 
and their associated measurement noise as a function of the survey parameters 
such as density and depth and for several popular models of structure formation. 
In addition, we have studied how well these kind of surveys can distinguish 
between different power spectra and found that, for the same number of objects, 
cone surveys are as good or better than full-sky surveys in distinguishing 
between popular cosmological models. We hnd that a survey with 200 — 300 
galaxy peculiar velocities with distance errors of 15% in two cones with opening 
angle of ~ 10° could put signihcant constraints on the power spectrum on 
scales of 100 — 300h“^Mpc, where few other constraints exist. We believe that 
surveys of the kind we describe her could provide a valuable tool for the study 
of large scale structure on these scales and present a viable alternative to full 
sky surveys. 


Subject headings: cosmology: distance scales - cosmology: large scale structure of the 
universe - cosmology: observation - cosmology: theory - galaxies: kinematics and dynamics 
- galaxies: statistics 
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1. Introduction 

One of the most important goals of cosmology is to determine the power spectrum of 
initial density fluctuations in the Universe. The primary tools in this endeavor have been 
redshift surveys, which have been used to probe scales < 100 h“^Mpc, and microwave 
background anisotropy measurements, which tell us about the power spectrum on scales 
> 1000 h~^Mpc (we use h = ifo/(100 km/s/Mpc) where Ho is the Hubble constant). The 
gap between these scales, where we have relatively little knowledge of the power spectrum, 
is notable in several respects. First, theoretical considerations tell us that the power 
spectrum should have a maximum in this region at a scale corresponding to the horizon 
size at the time that the Universe became matter dominated. While redshift surveys have 
hinted that the power spectrum does indeed turn over at a scale ~ 200 /i“^Mpc {e.g. Fisher 
et al. 1993, Feldman et al. 1994), the errors inherent in the measurement are too large to 
be dehnitive. Further, several recent studies have suggested that there might be much more 
power on these scales than is usually assumed in models of large scale structure formation 
(Broadhurst et al. 1990, handy et al. 1996, Ainasto et al. 1996.) If these suggestions are 
correct it will have important ramihcations for our understanding of how structure formed 
in the Universe. 

There are several observations planned or in progress that will attempt to study the 
power spectrum in this regime. From above, there are the small angle microwave anisotropy 
measurements. From below, the next generation redshift observations should signihcantly 
extend the range over which surveys can reliably determine the power spectrum. However, 
it may be a decade or longer before these ambitious and complex projects will produce 
measurements of the power spectrum on scales 200 h ^Mpc. 

A third method of measuring the power spectrum on large scales involves the study of 
large scale motions of galaxies. There were several recent attempts to determine the power 
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spectrum from observations of peculiar velocities. Kolatt & Dekel (1997) use the POTENT 
reconstruction of the density field to determine P(k). Zaroubi et al. (1996) estimated the 
power spectrum from the Mark III catalog of peculiar velocities using Bayesian statistics. 
One advantage of this method is that the large scale velocity held probes the matter 
distribution in the Universe directly, and not merely the light distribution as redshift 
surveys do. However, the errors in velocity estimates are typically some fraction of the 
distance of the sample points, which in the case of distant objects can mean that the 
errors are larger than the peculiar velocity being measured. This is partially rectihed by 
measuring only the lowest moment of the velocity held, namely the bulk how. Since the 
bulk how is, in a sense, an average of the velocities in the sample, its error is reduced over 
that of an individual measurement by the square root of the number of objects. 

Two recent ehorts to measure the bulk velocities of large volumes have resulted in 
contradictory conclusions. Lauer & Postman (LP 1994) found that the inertial frame 
dehned by Abell clusters within 15,000km/s exhibits a bulk velocity of ~ 700 km/s with 
respect to the cosmic microwave background (CMB) rest frame. More recently, Riess et 
al. (RPK 1995) used type la supernovae to measure the bulk velocity of a similar volume; 
their results are consistent with their sample being at rest relative to the CMB rest frame. 
In recent analyses (Feldman & Watkins 1994, Strauss et al. 1995, Jahe & Kaiser 1995, 
Watkins & Feldman 1995) it was shown that both power spectra from structure formation 
models as well as results from redshift surveys are inconsistent with the LP measurement 
at the 2 — 3(7 level, whereas they are quite consistent with the RPK result. Furthermore, 
the RPK and LP results seem to be inconsistent with each other at a high conhdence level. 

Abell clusters and type la supernovae are rare, thus the LP and RPK samples are 
by necessity quite sparse. This limits their ability to accurately measure moments of 
the velocity held beyond the bulk how. In addition, obscuration by our Galaxy strongly 
constrains how accurately these surveys can ascertain components of the bulk how in the 
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plane of the Galaxy. Thus the only component of the bulk flow that LP and RPK are able 
to report with a reasonable signihcance is that along the Galactic poles. In this Letter 
we propose an alternative approach to gathering full sky surveys of peculiar velocities; we 
explore the possibility of probing the power spectrum on large scales by measuring the 
velocities of galaxies in only small patches of the sky. While this survey will necessarily 
be sensitive to only one component of the bulk flow, the increase in the density of objects 
will allow more accurate measurement of higher moments of the velocity held. An added 
advantage of this approach is that the higher moments probe a different range of scales 
than the bulk how, making it possible to constrain the power spectrum more precisely. 

We examine several factors necessary for the design of a useful peculiar velocity survey. 
We concentrate on surveys which cover two helds directly opposite to each other. These 
surveys can probe larger scales and are less susceptible to radial biases than one-sided 
surveys. We construct mock catalogs for these surveys and show how their sensitivity 
depends on depth and the number of objects measured. In particular, we calculate the 
expectation values and expected errors associated with the three most easily measured 
moments assuming several diherent power spectra. 


2. Analysis 

As stated above, individual velocity measurements are too noisy to allow us to 
accurately map the peculiar velocity held v{r). Instead we expand v{r) in a given region in 
terms of its moments (Kaiser 1991, Jahe & Kaiser 1995), 

Vi{r) =Ui + r-jPij + rjVkQijk + ■■■ ( 1 ) 

Here, Ui and pij are usually referred to as the bulk how and the shear tensor respectively. 
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In practice we can measure only the radial component of velocities. Thus if our objects 
lie in a cone around a direction given by R, we will be most sensitive to the component of 
the velocity along R, vji = v ■ R. For this situation it is sufficient to model the velocity as 
being entirely along the R direction and depending only on f ■ R, giving us a much simpler 
expansion, 

VR{f- R) = U, + U2{f- R/Ro) + U^{r- R/Ro)\ ( 2 ) 

where we have introduced the arbitrary scale Ro in order that the constants Ui will all have 
the same units. In most of the analysis below we choose Rq = 100/i“^Mpc. The constants 
Ui represent the one component each of Ui, pij, and qijk to which these surveys are most 
sensitive. 


In terms of this model, the estimated line-of-sight velocity Sn measured for 
the nth galaxy at a position rh can be written as Sn = Y^i=iFn,i Ui-\- en, where 
Fn,i = ifn ■ R){fn ' R/Roy~^. Here we assume that the noise is drawn from a Gaussian 
with zero mean and variance where is the estimated uncertainty in the 

measurement of the line-of-sight velocity and cr* is introduced to account for contributions 
to the velocity of the galaxies in the survey arising from nonlinear effects as well as from 
the components of the velocity held that we have neglected in our model (Kaiser 1988). We 
shall take a* = 400 km/s for all of our calculations, although in practice this choice makes 
little difference given the large values of Cn for most of the galaxies that we will consider. 


Given a sample of objects with positions r„ and line-of-sight velocities Sn, the maximum 
likelihood solutions for the constants Ui are given by 


U, = (A-') 


-G.. V 

n 


Fn,j Sn 


.O': 




( 3 ) 


where the matrix An = -rS—AA', here and below we use the Einstein summation 
convention for repeated indices. 
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The theoretical expectations for the constants Ui can be expressed in the form of a 
covariance matrix, Rij = (UiUj) = where is the contribution from the 

velocity held and is the noise. The matrix is obtained by convolving the power 

spectrum with a window function, 

d;’ = ^/ p^imUk) <ek = j p[k)wf,[k) dk, (4) 

where we have used the fact that in linear theory there is a simple relationship between the 
velocity and density power spectrum, P{v){k) = {H‘^/k‘^)p{Qo)P{k), with f{flo) ~ 



Once we have calculated Rij for a given power spectrum, we can construct a 
statistic = Ui Uj for the three degrees of freedom of the measured moments. This 
statistic can be used to assess the compatibility of the power spectrum with the measured 
Ui- We note that the this statistic properly accounts for the correlations between the 
moments, which can be important for a sparsely sampled survey. Alternatively, given 
a parameterization of the power spectrum, Rij can be used in a likelihood analysis to 
determine the most likely values of the parameters (see, e.g. , Jaffe & Kaiser 1995). 


In designing a survey, it is important to determine how well it can distinguish between 
different power spectra. We can quantify this by calculating the expectation value of the 
for power spectrum PA{k) given that the true power spectrum of the Universe is power 
spectrum Psik): 

(xab) = mBURXhi(UB),) = {RXk,mBUUBh) = (.RXkARBk, = TtIrARb), (7) 



where subscripts denote quantities calculated assuming power spectra Pa and Pb- The 
value of (xab) can be used to estimate the confidence level at which one can rule out power 
spectrum Pa in a universe with power spectrum Pb given a particular catalog of objects 
and assumed velocity errors. 


We constructed mock catalogs of hxed number of galaxies restricted to two cones 
on opposite sides of the sky. We found that the results were roughly independent of the 
opening angle a of the cone as long as a < 30°, giving some flexibility in the design of a 
survey. The results we report in this Letter are for a = 10°. Galaxies were distributed 
radially using a selection function obtained by fitting an analytical function (cf. Feldman et 
al. 1994) to the radial distributions of various magnitude limited Tully-Fisher surveys (we 
used different surveys from the Mark III Catalog of Peculiar Velocities [Willick et al. , 1995, 
1996]). We adjusted the depth of the survey by varying the location of the peak of the 
distribution function, Ro, which is directly related to the magnitude limit. For dehniteness, 
we assumed the error in the measured velocity of a galaxy to be 15% of its distance, an 
error commonly reported for Tully-Fisher distance measurements. 


Given a mock survey, we can calculate the window functions for the three moments 
and determine which scales they are sensitive to. Assuming a power spectrum, we can also 
calculate the expected values of Ui and their associated noise. It is useful to dehne the 
parameters 


Qi 


d('^) 


( 8 ) 


that indicate how accurately a given survey can measure each of the three moments. Ideally 
we would like Qi > 1 for each of the three moments. The geometry and density of the 
survey can be optimized to meet this requirement with the minimal observational effort. 
Although Ui can in principle be correlated, we found that these correlations are small. 
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In Figure 1 we see a series of normalized window functions from three mock surveys 
with different values of Ro, each with the same number of galaxies (n = 200). While the 
window function for the bulk flow (f/i) has a maximum at fc = 0 , the window functions 
for U 2 and t /3 are peaked at smaller scales. Thus we see that the higher moments probe a 
different region of the power spectrum than the bulk flow. Since we are looking only at a 
single component of the velocity along the direction of the survey, the angular distribution 
of the objects is relatively unimportant and edge effects are negligible. 

In the bottom panel of Figure 1 we show the power spectra we used to calculate 
expectation values for the moments; standard COBE normalized CDM (Bardeen et al. 
1986) and HCDM (Klypin et al. 1993) as well an = 1 CDM normalized the observed 
abundance of clusters (see, e.g. Eke et al. 1996) and ACDM normalized to both cluster 
abundance and COBE. 

In Figure 2 we show the Qi parameters (Eq. |[) for the three moments for a survey with 
Ro = 100h“^Mpc (effective depth 200h“^Mpc) using the four power spectra described 
above. We see that we need some 200 — 300 galaxies in our survey to get Q* ~ 1 for the 
power spectra we consider that produce the largest velocities, i.e. COBE normalized CDM 
and HCDM. 

In the upper panel of Figure 3 we show the three degrees of freedom for 

distinguishing between HCDM and ACDM in a 10° cone survey with Ro = 100h“^Mpc as a 
function of the number of objects. The upper curve is for HCDM assuming that the true 
power spectrum is ACDM; here the > 3 since we are testing a model with excess power 
relative to the actual spectrum. The lower curve is for ACDM assuming the true spectrum 
is HCDM, which will give < 3. From these x^ values, one can assign a conhdence level 
at which the models are ruled out. In general, one can rule out spectra with excess power 
at a higher conhdence level. For comparison, we show the same quantities for a full-sky 
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survey (with the zone-of-avoidance removed) using the three bulk flow components as the 
three degrees of freedom; these surveys have the same radial distribution, number of objects 
and errors as the cone surveys. In the lower panel we show the for COBE normalized 
CDM given cluster normalized CDM and vice versa. From the figure we see that the three 
moments calculated for the cone survey do better than the three bulk-flow components of 
the full-sky survey at distinguishing between the models for the same number of observed 
galaxies. 

A more complete comparison between full-sky and cone surveys would include 
information from higher moments of the full-sky survey. However, for the surveys and 
power spectra we have considered, the signal to noise of these higher moments is small. 
Results from analyses including the highest signal to noise moments suggest that cone 
surveys continue to be as good or better at distinguishing spectra. 


3. Conclusions 

In this Letter we have explored the properties of proper distance surveys that cover 
small fields in two opposing directions. Our analysis exploits the fact that a small area 
survey can measure some of the moments of the velocity held much more accurately than 
a full sky survey with the same number of objects. We have shown how to expand the 
velocity held in moments and constructed a test useful for constraining models. In order 
to get a “signal to noise” of unity for the three lowest moments, we found that a survey 
of ~ 200 — 300 galaxies is needed if we assume distance indicators accurate to about 15% 
of the distance in a cone survey with opening angle of (9(10°) and consider survey depths 
~ 200/i“^Mpc. We have also shown that cone surveys are as good or better than full-sky 
surveys at distinguishing between cosmological models. These surveys could put important 
constraints on the power spectrum on large scales with only a modest observational effort. 
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and thus could provide a valuable tool in probing scales that have been up to now largely 
beyond our scope. 
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Fig. 1.— In the top three panels we plot the window functions associated with the moments 
f/i, f /25 and f /3 for different values of Ro, the location of the peak of the radial selection 
function for mock surveys of 200 galaxies. In the bottom panel we see the power spectra we 
used in our analysis. To find the expectation values of the different moments we integrate 
the product of the window functions of the top three panels with each of the power spectra. 
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Fig. 2.— Here we see the values of Qi for different power spectra as a function of the number 
of objects in the survey. Qi rises with the number of points as expected as the noise inherent 
in the measurement becomes smaller with increased density. 
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Fig. 3.— The fo^' 3 degrees of freedom for HCDM normalized to COBE given a Universe 
with ACDM, and vice versa, for 10° cone surveys, with the three lowest moments and for a 
full-sky survey (with the zone-of-avoidance removed) using the three bulk flow components 
as the three degrees of freedom, as a function of the number of objects. The full-sky surveys 
have the same radial distribution and errors as the cone surveys. We also show the 
COBE normalized CDM and cluster normalized CDM. From the hgure we see that the three 
moments calculated for the cone survey do better than the three bulk-flow components of 
the full-sky survey at distinguishing between the models for the same number of observed 
galaxies. 
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